Controllable tt junction in a Josephson quantum-dot device with molecular spin. 
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I. INTRODUCTION 



■ We consider a model for a single molecule with a large frozen spin sandwiched in between two 
BCS superconductors at equilibrium, and show that this system has a tt junction behavior at low 
temperature. The tt shift can be reversed by varying the other parameters of the system, e.g., 
temperature or the position of the quantum dot level, implying a controllable tt junction with novel 
application as a Josephson current switch. We show that the mechanism leading to the n shift can 
be explained simply in terms of the contributions of the Andreev bound states and of the continuum 
of states above the superconducting gap. The free energy for certain configuration of parameters 
shows a bistable nature, which is a necessary pre-condition for achievement of a qubit. 
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^ • Molecular spintronics is a promising new domain, at the convergence of two challenging disciplines. On the one 

■ hand there is molecular electronics, where single molecules are used to create electronics devices at the nanometric 
. scale with unique properties, while on the other hand there is spintronics, where the spin of the electron is used as the 

relevant quantity in place of the electronic charge. The latter allows us to take advantage of the unusual properties 
'"^ , of spin, like a long coherence time. It is in this context that we consider in our work the equilibrium properties 
^ ' of a molecule with a large magnetic moment placed between two superconductors, when a Josephson current flows 
O ' between the two superconductors through the molecule. The focus will be on the effect, on the Josephson current, 
of the coupling between the spin of the electrons producing the current and the molecular spin. We will show in 
particular that when this spin coupling is large enough, the superconducting junction behaves as a tt junction, with 
. a reversal of the Josephson current compared to the case without spin coupling. 
^ ■ Of great importance for molecular spintronics are molecules possessing a large spin, or "single molecule magnets" . 
9^ Such molecules can now be synthesized, for example the molecule Mnl2ac, which has a ground state with a large 
spin S = 10, and a very slow relaxation of magnetization at low temperature.^ This slow relaxation is due to a high 
\ anisotropy barrier for the magnetization, around 5.6 meV for Mnl2ac. For the system we are considering, this is a very 

■ large energy, as the typical energies in our system (temperature, coupling to the electrodes, etc.) are at most of the 
\^ ' order of the superconducting gap (which is 0.1 meV in Aluminum for example). This motivates our choice to take the 

, molecular spin as a fixed quantity, which will act as a local magnetic field for the electrons going through the molecule. 

■ Note that other systems, where the spin is not fixed, involving for example superconducting transport through fuUcrcne 
j3 . molecule doped with magnetic impurities have been considered experimentally and theoreticallyi^ Concerning the 

electronic transport across the molecule, we model the molecule as a single resonant level, i.e. a quantum dot. As we 
will be interested in the regime of good transparency between the molecules and the superconducting electrodes, we 
will neglect in this work the electronic interactions on the resonant level."* 
Q I The main result of our paper is to show that, when the coupling to the molecular spin is large enough, the system 
O ' shows a TT shift. A reversal of the super-current in a Josephson device and the free energy having a global minima 
at phase difference tt is referred to as tt shift and a Josephson junction displaying this is termed a tt junction.'^ The 
TT junction has potential applications in superconducting electronics, in quantum logic circuits as switches and are 
an integral part of superconducting phase qubits. We will also show that this tt shift can be controlled by the other 
parameters of the system (position of the dot level, temperature, coupling to the electrodes, etc.), allowing to reverse 
the TT shift and recover a standard Josephson junction. 

Generally in works related to tt junction behavior the bound state current, which is due to current carrying Andreev 
bound states formed between the two superconductors, is investigated while the contribution from the continuum of 
states above the superconducting gap is ignored^. There are good reasons for doing so, since the continuum contribution 
is generally much reduced compared to bound state current, especially in the limit of a long junction or a very short 
one. However recent works have shown that the continuum current cannot be ignored ^iSiiSiii^ especially in a Josephson 
junction which is neither very short nor very long. In this work, we calculate explicitly the contribution from the 
continuum, and we show that the in the presence of a large coupling to the molecular spin, the continuum current is 
essential to understand the tt junction behavior. In some regime, the bound state current can even vanish, and the 
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FIG. 1: Our model system 



continuum current is then the only contribution. We will show also that, with some fine-tuning of the parameters of 
the system, the system can be in a bi-stable state, where the (j) = and cj) = tt state are equally stable; this bi-stability 
is a necessary condition for a possible qubit implementation. 

The rest of the work is organized as follows. The next section deals with a short history of the 7r-shift as seen in 
Josephson junctions and with the possible applications which such behavior may have. Section III is devoted to the 
derivation of the Josephson current when coupling to the molecular spin is present. In section IV we use the formulae 
obtained in section III to show the behavior of the Josephson current as a function of the coupling to the molecular 
spin, and of the other parameters. We give a detailed explanation of the mechanism leading to the tt shift. In section 
V we discuss some potential applications of our system, first as a Josephson current switch, then as a superconducting 
qubit. Section VI is devoted to concluding remarks. 

II. BRIEF HISTORY OF vr-SHIFT 

In order to show how our work and results differ from existing works on tt junctions, we find useful to give 
a very ^ short history of the 7r-shift. tt junctions were first proposed theoretically by Bulaevskii and coworkers in 
Ref. 11]. They considered a tunnel junction with magnetic impurities in the barrier. In this system spin-flip tunneling 
leads to a formation of tt junction. They also predicted that a super-conducting ring containing a tt junction could 
generate a spontaneous current and a magnetic flux opening the way for experimental detection. Spin flip tunneling 
in superconductor-quantum dot-supcrconductor(S-QD-S) system has also been shown to give rise to a tt junction 
behavior as in Refs. [12..13tJ4] . It was Kulik who in 1966 was the first to discuss the spin-flip tunneling through an 
insulator with magnetic impuritiesj^^ The spin-flip tunneling is predicted to dominate the Josephson current when 
spin on the quantum dot is non-zero. In S-QD-S junction, changes in the sign of the critical current could be observed 
as a function of the quantum dot gate voltage which controls the occupancy of a quantum dot. Due to this gating 
capability one has more control over the magnetic state of a barrier in S-QD-S junction compared to a magnetically 
doped Superconductor-Insulator-Superconductor junction.— Superconductor- Ferromagnet-Superconductor(SFS) have 
also been shown to give rise to a tt junction behavior both in theoryiiii^ii^ as well as in experiments j ^° ^ The study of 
the superconducting tt state sheds more light on the coexistence of superconductivity and ferromagnetism in general 
and is also important for superconducting electronics.^^ Generally with increase in the strength of the exchange field 
the TT shift is observed with a reversal and suppression of the super-current. In SFSFS systems with the ferromagnets 
in anti-parallel alignment, however, with increase in the strength of the exchange field an increase in super-current is 
observed.— In SIFIS and SFcFS, where c denotes a constriction, structures also such tt junctions have been observed?^ 
Recently triplet superconductor-ferromagnet-triplet superconductor junctions have been predicted to have potential 
applications as current switches.-^"* In contrast to SFS systems, tt junction behavior in SNS systems occurs due to the 
creation of a non-equilibrium distribution of electrons in the barrier via a control channel. Thus in these systems the 
TT-junction can be controlled via a voltage applied to the control channel this makes such devices ideal for them to be 
used in superconducting digital circuits, especially as a phase inverter, i. e., tt-SQUID^^ for complementary Josephson 
digital devices. Further application of tt junctions as candidates for engineering quantum bits have been predicted;^ 
Finally, tt junctions have also been theoretically predicted and experimentally observed in superconducting d-wave 
j unctions 
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III. DERIVATION OF THE JOSEPHSON CURRENT 
A. Model Hamiltonian 

The Josephson Current(/,7) can be calculated from the derivative of the free energy(F) with respect to the phase 
difference ((/)) across the superconducting leads Ij — 2^, in equilibrium. The free energy in turn is defined as 
F — —kTlnZ, where Z is the partition function of our system. Thus calculating Z is the first step in calculating the 
Josephson current in our system. The full Hamiltonian of our system is written as below: 

H^Hd-s+ J2 Hj+Ht, (1) 

j=L,R 

where H^^s defines the Hamiltonian of the quantum-dot molecule with spin, Hj represents the superconducting 
leads, while Ht denotes the tunneling part. The dot-spin Hamiltonian is: 

HD-s^eY,dU<r + JS{d\d^-d\di) (2) 

(T 

where d„,d\ are the electronic operators in the dot, e is the energy level of the dot, J the coupling between the 
molecular spin S and the electronic spin on the dot level. The coupling term comes from the exchange interaction 
JS • s, where s is the electronic spin on the dot, but as explained in the introduction, the molecular spin is fixed in 
our system, and we chose the spin quantization axis along the spin orientation. In the superconducting Hamiltonian 
it is convenient to perform a gauge transformation which removes the phase from the order parameteri^^ Thus, 

J, \ ^i(-fc),i / 

and finally the tunnelling part can be written in the standard form with a hopping parameter tj determining the 
transfer properties of the junction. The Pauli matrices mentioned in the above equation are matrices in particle-hole 
space. The effect of the gauge transformation on the tunnel Hamiltonian is the appearance of a phase dependence in 
the hopping parameter 

HT = Y.^],T,d + h.c., (4) 

with Tl = tLaze'^"''^/^,Tii — i^CTze"*'^^"^/^, where (j) is the phase difference between the superconducting leads, and 
tj is the tunnelling amplitude between the jth lead and dot. 

B. Neglecting Coulomb interaction 

In the present work we chose to neglect the charging energy of the molecule. Molecular electronics transport 
calculations are typically concerned with two limits: either the limit where the tunnelling rate from the molecule to 
the leads dominates over the charging energy, or the opposite limit of strong Coulomb blockade which can be dealt 
in the incoherent regime or in the coherent regime. The validity of each regime depends on the transparency of the 
tunnel barriers connected to the molecule, and on the capacitances seen by the molecular quantum dot with respect 
to the leads. So far in the literature, many efforts have focused on the Landauer-Buttiker description of molecular 
electronics transport In some instances^ this approach is supplemented by taking Coulomb interactions in an 
effective manner using density functional theory, but the Green function which is used to compute the transmission 
coefficient is determined from a specific electron configuration, as for an effective one-electron Coulomb potential. 
Here we adopt a Hamiltonian approach to compute the current, but the basic assumptions for neglecting Coulomb 
blockade effects also apply. This is justified as follows. 

In fact this stems from the fact that the charging energy of the molecular quantum dot is assumed to be small 
compared to the the escape rate of the electrons from the dots to the leads, later to be referred to as F. This regime 
triggers a substantial broadening of the dot level, as for instance was discussed in Ref.flU where a molecule was 
attached to metallic substrate, being effectively "metalized" . Qualitatively speaking, if the transparency from the dot 
to the leads is close to unity, the time scale which characterizes the lifetime of electrons on the molecular quantum 
dot is so short that Coulomb blockade effects do not have time to operate. 
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These qualitative arguments are substantiated by theoretical works on the Coulomb blockade in the presence of 
highly transmissive barriers, which were carried out more than a decade ago. Ref.jssj studies the behavior of a 
quantum dot embedded between two point contacts, and finds a crossover to a regime where charge fluctuations in 
the dot are dominant, therefore wiping out charge quantization (Coulomb blockade) effects. The temperature which 
characterizes this crossover is of course proportional to the dot charging energy, but it also goes to zero in the limit of 
ideal transmission. Ref. 34] considers a quantum dot connected to single barriers, and shows that the energy of the dot 
undergoes (Coulomb blockade) oscillations as a function of gate voltage as long as the transmission coefficient of the 
barrier, which isolates the dot, does not approach unity. On another note, Ref. [3^ uses a path integral framework^ 
to describe a quantum dot with arbitrary barriers. Previous results^ concerning charge fluctuations are recovered, 
but more importantly it is shown that the effective charging energy is exponentially reduced at ideal transmission. 
This latter result applies, granted to a dot coupled to several channels, but these features are expected to survive for 
a a dot coupled to a single, highly transmissive channel. 

Note that in the above, it is sufficient for only one of the two contacts to have a large capacitance in order to be 
able to neglect Coulomb effects. The present point of view is consistent with recent works on molecular electronics 
issues where phonons are involved^, but electron-electron interactions are neglected nevertheless. 

Finally, we stress the fact that there exist actual experiments in molecular electronics, which can achieve the 
high transmission conditions which are assumed in the present work. It has been demonstrated that break junction 
geometries^ can achieve close to ideal transmission: for a Hydrogen molecule which is sandwiched between Platinum 
electrodes, one observes a conductance quantization plateau at 2e^ /h, when bringing the two electrodes together, as 
expected from a monovalent metal as Hydrogen (in contrast, a pure platinum junction yields steps at 4e^//i). This is 
the direct evidence of single, perfectly transmitting channel. To summarize, in these works, which have been extended 
to study the effects of vibrations on the molecule^, no Coulomb blockade effects show up at all. For the present 
study, we thus assume that our molecule is placed under the same conditions as in these break junction experiments. 
Note that high transmission conditions have also been obtained with Carbon nanotubes^i^. 



C. Effective action 



To calculate the partition function we use the path integral approach. In this method the partition function is given 
by: 

Z = JlliV^jkV^jkVdVd)e-^^ (5) 

jk 

Z is written as a functional integral over grassmann fields for the electronic degrees of freedom (^P, 5*). The Euclidean 
action reads: 



Sa^Sd+ / dT[V ^jk{T)(dr + Ck<J, + ACT,)*jfc(r) + Ht{t)] 
Jo ,,, 



/3 

lr[] 

Jk 

/3. 



P is the inverse temperature, and Ht{t) ~ Yljk^ jk{T)Tjd{T) + h.c. while Sd = Jq dT[d{dr + euz + JS)d]. After 
integrating out the leads we get 



Z = jvdVde-^^" with Sejj =Sd- j^drdr' d[T)t{T-T')d[T') (6) 

where t{r) = Y.j=L,RT}G{r)Tj and G(t) = Efc(9r + S,kcr. + AfT,)-i,5(r). 

We perform a Fourier transform on the Matsubara frequencies (with Wn = (2n+ l)7r//3): S{t) ~ i e^™""^ and 
G(r) = i X)«;„ e~™"'^G^,^, which gives for the Green function G: 

Gw„ = Idi y{S,){-iwn + ik<Jz + Act^)-i ~ + ^^^) ('') 

J ^wi + 

In the above equation, — X^fc ^('? ~ ^k) is approximated as a constant i'(O), the density of states at the Fermi 
level in the normal leads. This gives for the self-energy: 

= ^^'^ [iwn - A cos((/./2)a, - 7A sixv{(t>/2)ay] (8) 

a/a-^ — {iwy 
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with 7 



1 

VP 



ri+r„ 

— iWnT^ 



, r = + ^R,^L/R = ^■^^(O)^/^//^- We get finally for the effective action (introducing (i(r) 

.J 



^■A^u;„rfti)„ with = —iwn + ecr^ + JS* — 



(9) 



D. 



Andreev Levels 



The dispersion equation for the Andreev levels is given by the eigenvalues of the effective action in Eq. ([9|) (with 

iw = z) 



det 



z — euz — JS 



r/2 



VA2 - z2 



(z — A cos(0/2)(T2: — 7A sin(0/2)o'j,) 







which gives (introducing the parameter s = JS): 
Tz Tz 



2VA2 - z2 



e- s){z + 



2VA2 - z2 



r2A2(cos2((/./2) + 72 sin^(0/2)) 
4(A2 - z2) 







(10) 



(11) 



While this cannot be solved analytically in general, there are two limiting regimes where one can get an analytical 
expression of the solutions, giving two Andreev levels. For simplicity, we choose here e = and 7 = 0. 
case 1: T > A : 



z — Ei,2 — A cos 



/ ± cos 0/2 \ , /2s 
Arccos — = + Arctan — 



(12) 



case 2: T, s < A 



El.; 



s ± - cos((/./2) 



(13) 



In the general case, for arbitrary F, e, s(= JS) and A, we calculate numerically the roots, by transforming the l.h.s. 
of Eq. pT|) into a 8th order polynomial in z to get rid of the square roots, and then calculating the roots of this 
polynomial. We find that only two of these roots correspond to roots of Eq. (fTTjl (see also ref. Q), and that these 
two roots are real and belong to [—A, A]. There are thus always two Andreev bound states, as in the zero spin case: 
the effect of the spin term is merely to move these two states, but it does not introduce new bound states. 

In Fig. [21 we plot the two Andreev bound state positions as a function of the phase difference for four values of 
spin, \s\ = 0, 1, 2, and 4, for large transparency of the contacts (F = 4.0) and very low temperature {(3 = 1000). A = 1 
is taken as the unit of energy in our system, as in the rest of this work. The right panel in Fig. [2] corresponds to s > 
while left panel is for s < 0. We see that when the absolute value of s is increased, the two Andreev levels are pushed 
towards +A or —A. 

It might seem surprising that the Andreev bound states are different for +s and — s. Indeed, our physical system 
is invariant under the interchange of spin up and spin down electron combined to the exchange s — s; as the 
superconductors are the same under the spin interchange, our physical system has thus to be invariant under the 
transformation s — s- —s. As will be shown in the next sections, the total Josephson current, which is a physically 
measurable quantity, is invariant under this transformation (s —s). However, the expression of this total current 
in terms of the Andreev bound state current and of the continuum current, depends on ths sign of s (and so do the 
Andreev bound states) 

The different Andreev bound states obtained for s and for —s are thus two different ways to represent the same 
physical situation. The fact that we obtain one of the two possibilities for a given external spin can be traced back 
to our initial choice of the spinors for the superconductors (Eq. ([3])) and for the dot (Eq. ([4])): had we chosen the 
spinors defined with opposite spin (for example, d) — (d|,(i|), to be compared with the definition of Eq. (|1])), we 
would have obtained for s < the Andreev bound states shown here for s > 0, and vice- versa. One could also use a 
combination of the two possibilities of spinors, in order to get a spectrum of Andreev bound states (and bound state 
current) which is independent of the sign of s; in this case, the spectrum is composed of four Andreev bound states, 
which are precisely the two we obtained for s > plus the two for s < 0. In this paper, we have chosen to keep the 
spinors as defined in Eq. ^ and Eq. This choice will give us a particularly simple picture for the mechanism 
leading to the pi-shift (see below). 
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FIG. 2: The two Andreev bound states as function of phase difference. For s > in the right panel and for s < in 
the left panel. The labelling of the curves is as follows: spin s=0(black, solid line), 1.0 (red, dotted line), 2.0 
(blue, dashed line) and 4.0 (green, dot-dashed line). The other parameters are F = 4.0, 7 = 0.0, /3 = 1000, e — 0.0, 

and A = 1.0. 



E. Josephson Current 

The partition function after integrating out the {d, d} variables is given by- 

Z = jVdVde-^^" ='[[detMru,, (14) 
where A4nj„ is given in Eq. ([5]). The Josephson current then reduces to: 

1^ F^A^(l-7^)sin(0) _ 1 ^ 

/3 2^4[detA^„J(A2-(zu;„)2) - p 2^ ^'^^ 



Wn — — 00 



where the last equality defines the function /. 
Further, the Free energy is given by: 

i^ = -i^ln(detA^^J. (16) 

In the above equations, detAiw„ is the same as the l.h.s of Eq. (fTTI) . with iwn replacing z. 

^From the above equation, one can calculate the total Josephson current by summing over the Matsubara frequen- 
cies. However, we can transform the above equation in order to separate explicitly the contributions of the Andreev 
bound states and of the continuum, which are physically meaningful. In order to calculate these contributions, we 
take advantage of the fact that the Matsubara frequencies are the poles of the Fermi function nf{z) [4l|. We then 
consider the integral I = j^f{z)nf{z), where the function f{z) is defined in Eq. p5| . The function f{z) as seen 
earlier has two poles on the real axis between —A and A (these are simply the two Andreev bound states, for which 
det — 0). Further, because of the square roots terms in the det A^, it has branch points at z = ±A; we have 
chosen to place branch cuts on the real axis, for z G [A, 00] and z G [—00, —A]. We thus chose the contour C as two 
large semi-circles plus parts going around the branch cuts. We illustrate the contour, poles and branch cuts in Fig. 
3. Thus integral I can be broken into the sum of the contributions from the large circle D of radius R, the two small 
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FIG. 3: The contour includes poles from the Fermi function, two Andreev bound states and the contribution from 

the branch cuts. 



circles at ±A, denoted by di^2 of radius e, and the contribution from the branch cuts. Therefore, 

1=77-1 f{z)nf{z)dz = 77- I dz f{z)nf{z) + ^ [ dz f{z)nf{z) + ^ [ dz f{z)nf{z) 

1 /"^ 

1 

+ ir^l™/ [/(z + ze) - /(z - ie)]n/(z) 



(17) 



The integrals over D and c?i_2 tend to zero as i? — > oo and e ^ 0. The last two terms in the above equation define the 
contribution from the continuum to the current, which we denote as /c. From Cauchy's residue theorem the integral 
I can also be evaluated as follows: 

(a) Function n/(z) 

-^Y.w„ f'y^'^n), which is precisely /j (Eq. dTS])). 

(b) j[z) has 2 poles of its own, written E\ and E^- These gives the Andreev bound states contribution, which we 
denote by /(,. 

We have thus 



j^q-j- has poles at z = iwn-, with residue —1/(3. The contributions from these poles is thus: 



f I, J =h+Ic 

h - -nF{Ei)ves{f{Ei)) - np{E2)ves{f{E2)) 



Ic — lim — 


/>oo 


.-A 




/ dz + 


/ dz 


'^{f{z + ie)nf{z + ie)) 




J A 





(18) 



where res denotes the residue of the quantity in square brackets, and 3 stands for the imaginary part. We see from 
the equation giving If, that the contribution from each Andreev bound state is simply proportional to the occupation 
number npiEi) (z = 1, 2) of this level. Equation (fT8|) is the central result of this work, which we have used to calculate 
numerically the lb and Ic curves shown in the following. 



IV. RESULTS 



In this section, we show the results obtained numerically for the Josephson current using equation (|18p . We will 
see that, when the coupling to the spin is strong enough, the junction behaves as a tt junction, and that the other 
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FIG. 4: The total Josephson current (/,7)as function of the phase difference for increasing spin. The tt junction 
behavior is clearly seen, (a) F = 0.2. (b)F — 2. (c)F = 10. The other parameters are: A = 1.0,(3 = 1000,7 = 0: ^-nd 

e = 0. 



parameters of the system allow a control of this tt junction. We will also show that the mechanism leading to this 
TT junction behavior can be understood simply in terms of the bound states and continuum contributions to the 
Josephson current. 



A. Total Josephson current and free energy: The tt shift 



In Fig. m we plot the Josephson currcnt(/j) as a function of the phase (p, for different values of the spin coupling 
s = JS. We clearly see a 7r-junction type behavior as the magnitude of the spin coupling is increased. One clear 
inference from Fig. 3] is that the strength of the spin coupling required to engineer a tt junction behavior increases 
with the interface transparency F. We get another point of view of this vr shift in Fig.[5l where we plot the free energy 

as a function of the phase 0. We see that with increase in spin coupling strength the transition from the to the 
TT phase is clearly marked: the absolute minimum of F shifts from = to = tt. 

This figure also brings out other features, namely the 0' and tt' phases. As is evident from Fig. 5, the labelling 
of the respective junctions as 0,0', tt' and tt configurations, follows from the respective stability of = and </> = tt 
configurations. For a 0(7r) junction, only — 0(0 = tt) is a minimum of F{(j>). For the other two cases, both = 0, tt 
are local minima and depending on whether </> = 0(0 = tt) is the global minimum, one has a 0'(7r') junction. Of 
particular interest is the bistable junction, in which both = and = tt are global minima; note that this bi-stability 
is a necessary precondition for the realization of a Josephson junction qubiti^ 

The total Josephson current and the free energy are invariant with respect to the change of sign of the spin (s —f —s). 
As explained in part B of the previous section, this is to be expected from the invariance of the system under the spin 
up - spin down exchange. Technically, it can be seen on Eqs. (|15p - (|16p . using Lu_(^n+i} — ~^n- Note that, as explained 
before, the Andreev bound states are not invariant under s — > — s, and thus the Andreev bound states current and 
the continuum current are also not invariant. 



B. Mechanism of the vr-shift 



The ability to distinguish, in the Josephson current, between the contributions from each Andreev bound state 
and from the continuum (see Eq. (fTS]) ) provides us with a simple picture for the mechanism leading to the tt shift 
for large spin coupling. Note that the picture we obtain depends on the initial choice of spinors (see the discussion 
at section III D); the choice we make here allows us to get a very simple picture. In a few words, the effect of the 
spin coupling is to reduce or suppress the Andreev bound states contribution, and to give more importance to the 
continuum contribution, and this leads to the tt shift. With more details, the effect of the spin coupling on the bound 
states current can be understood from Eq. and Fig. [21 For s = 0, we see in Fig. [2]that there is always one bound 
state below the Fermi level, and the other one above. As the contribution of a bound state to the Josephson current is 
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FIG. 5: The Free energy (F) as function of the phase difference (<i>) for increasing spin from top to bottom. The tt 
junction behavior is clearly seen. The other parameters are: A = 1.0, F = 4.0, 7 = 0, /3 = 1000 and e = 0. 



proportional to the occupation number np{E.i) of this bound state (Eq. (|18[) ). this means that we have only one bound 
state contributing to the current, and this contribution appears to be much larger than the continuum contribution. 
With a large positive spin coupling s, we see on Fig. [5] that both Andreev bound states are above the Fermi level, 
which means their contribution to the Josephson current vanishes; while for a large negative spin coupling s, we see 
that both bound states are below the Fermi energy, which means they both contribute to the Josephson current, and 
this reduces the total bound state contribution as the respective contributions of the two bound states have opposite 
signs. Note that the total Josephson current is independent of the spin coupling sign, but that for large s > there is 
only the contribution from the continuum, while for large s < there is a combinations of the bound states and the 
continuum contributions. This explanation is illustrated on Fig. [SI where the contributions of the bound states and 
of the continuum are plotted for different values of s. 

The origin of the continuum current -which is non-zero even at zero temperature- is due to the phase difference 
between the two superconductors, which breaks the symmetry between the left and right-moving quasi-particles42ii^ 
One can draw an analogy with persistent currents flowing in normal metal rings at zero temperature. In normal metal 
rings the flux breaks the symmetry between clockwise and ant i- clockwise moving electrons inducing the persistent 
current. At zero temperature all states below the Fermi energy are filled, still then the persistent current is non-zero4i 

We also observe that the continuum current generally flows opposite to the bound state current. This observation 
is in agreement with that of other worksi^ 

We finally add that one can also understand the fact that the full current is the same for +s and — s but with 
different contributions from bound states and continuum by using electron hole symmetry. At first sight, electron 
hole symmetry only holds when the dot level coincides with the superconducting chemical potential. We first discuss 
this case, and then we address below the case where a gate voltage shifts the dot level away from this location. 

Consider the case of negative coupling (s < 0) in Fig. [2^. From the electron point of view, occupied states below 
the Fermi level have a continuum contribution to the current and a bound state contribution. From the point of view 
of holes, which occupy all states above the Fermi level, there is a continuum contribution to the current of occupied 
hole states while the contribution of bound states above the Fermi level diminishes with increasing \s\. For sufficiently 
large |s|, the two bound states are below the Fermi level and the holes cease to have a bound state contribution. On 
the opposite, for positive coupling in Fig. ^Bp, the role of electrons and holes is reversed: from the electron point of 
view, the bound state contribution is reduced -and eventually vanishes - when increasing s > 0; from the hole point 
of view, holes occupying both the continuum and the Andreev bound states above the Fermi level. We thus see that 
for positive (negative) s, the role of electrons and holes is reversed, and the electron hole symmetry can explain why 
physically observable quantities are invariant under the substitution s — *■ — s. 
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FIG. 6: The Andreev bound state 
current (/;,), continuum current (/c) and 
total Josephson current (/j) as function 
of phase difference (j). The individual 
Andreev bound state contributions 
(-^hi, Ib2) are plotted in the insets. 
Central figure is for spin s — 0.0. The 
two figures to the bottom of the central 
figure depict the case of s < 0, the first 
is for spin s = —1.5 and last one on the 
bottom is for s — —2.1. On the top of 
the central figure the two figures depict 
the case of s > 0. The individual figures 
plotted are for spin s = 1.5 and lastly 
for s = 2.1. The other parameters are: 
A = 1.0, r = 4.0, 7 = 0.0, P = 1000 and 
e = 0. 
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FIG. 7: (A) The Andreev bound state current (/;,, dotted line), the continuum contribution {Ic, sohd hne) and the 
total Josephson current (/j, dashed line) as function of dot level e, for positive coupling. The n junction behavior is 
clearly seen as dot level is varied, in the insets the individual bound state contributions (/f,i, /f,2) are plotted for s < 
and positive s > (for s < the continuum and total bound state currents are also plotted). The total Josephson 
current (//) is identical for s < and s > 0. The other parameters are: A — 1.0, s — 2.1, ip — tt/2, T — 4.0, /3 — 1000 
and 7 = 0. (B) The Free energy as function of the phase difference (0) for increasing dot energies, and (C) Andreev 
levels for both negative as well as positive coupling. The other parameters for B and C are: 
A = 1.0, s = 2.1, r = 4.0, (3 = 1000 and 7 = 0. 



Next consider the case where the dot level does not coincide with the chemical potential (e ^ 0). In any normal 
metal devices, this indeed breaks electron hole symmetry, nevertheless we wish to point out that because we are dealing 
with a superconducting system the above argumentation still holds. When dealing with an arbitrary superconductor- 
normal metal- superconductor (SNS) junction, Andreev bound states are understood from the conversion of electrons 
into holes and vice versa at the superconducting junction. An alternative picture is to say that two electrons, one 
above the chemical potential with energy E, and one below with energy —E are transfered through the normal region 
from one superconductor to the other. For our quantum dot setup, calculations of the Josephson current could for 
instance be performed using T-matrix formalism (to all orders for an exact result), where information about energies 
appears only in the form of differences of energy levels [E — e for the positive energy electron and £^ -I- e for the negative 
energy electron). These two energy differences are invariant under the change e — > — e. A dot with say, e > will thus 
have the same Andreev bound states spectrum and Josephson current, as a dot with e < 0. This is indeed explicit in 
Eq. pip , which is invariant under the transformation e — e. In the presence of the impurity spin the bare dot level 
is split by a Zeeman like coupling (the exchange term), nevertheless the above statement (e — e symmetry in the 
energy differences) still applies: the Andreev spectrum is reversed between +s and — s, but the Josephson current is 
shown to be the same because it can be computed from the point of view of electrons (say, for -\-s) or, alternatively, 
of holes (for ~s). 

In conclusion, because the Andreev bound state is the same for positive and negative e, the symmetry -\-s — )■ —s 
also holds for physically observable quantities such as the current. The argument for e = describing the contribution 
of electrons and holes populations on the continuum and discrete levels is unaffected by the modification e 7^ 0. This 
fact is demonstrated in Fig. [71:, where the bound state spectrum is plotted as a function of flux for several values of 
e, for positive and negative s. 



C. Controlling the 11 shift 



A remarkable feature of our system is that the tt shift behavior can be controlled and reversed using the different 
parameters of the system. This is important for potential experimental implementations, as some of these parameters 
can be accessed relatively easily (one can for example move the dot level by using a gate voltage'*^), while the spin 
coupling is a fixed quantity which depends in the molecule used. Our results show that, when the spin coupling is 
large enough to have a tt junction, a change in any of the parameter of the system (dot level e, coupling to the leads F, 
asymmetry of this coupling 7, and even the temperature) makes it possible to have the system behave as a standard 
junction (going through any intermediate situation between tt and junction). Schematically, the mechanism for 
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FIG. 8: (A) Andreev bound state current {h), the continuum contribution (/c) and the total Josephson current 
(/j) as function of inverse temperature /3 for positive coupling. The tt junction behavior is clearly seen as inverse 

temperature f3 is varied. In the insets the individual bound state currents for the case of s > and s < are 
plotted. For s < 0, the continuum along with the total bound state current is also plotted. The total Josephson 
current Ij is identical for s > and s < 0. The other parameters are: A = 1.0, s = 2.1, e = 0.0, 7 = 0.0, = tt/2 and 
r = 4.0.(B)The Free energy as function of the phase difference at different temperatures. The junction is bistable at 
the crossover between — tt junction behaviors. The other parameters are: A — 1.0, s — 2.1, e = 0.0, 7 = 0.0, and 

r = 4.0. 



this can be understood along the same lines as the explanation given above for the tt shift: starting from a tt junction 
situation, where both Andreev levels are (for example) above the Fermi energy and thus do not contribute to the 
Josephson current, changing a parameter of the system can move the Andreev level positions, and as soon as one 
of the Andreev level goes below the Fermi energy, it gives an important bound state contribution which brings the 
system back to a junction behavior. This is illustrated on Fig. [71 where the dependence of the currents (panel A), 
of the free energy (panel B) and of the Andreev levels (panel C) as a function of the dot level position e. Similar plots 
are obtained when looking at the F or 7 dependence (not shown). 

The picture is a bit different when the temperature is changed, as there the Andreev levels do not move, but the 
Fermi functions become broader as temperature is increased, leading to a partial revival of the bound state current. 
This is shown on Fig. [8l where the dependence of the currents (panel A) and of the free energy (panel B) is shown 
as a function of /? = l/{kBT). Starting from low temperature (high (3), with a tt junction behavior (the total current 
/j is < 0, and the free energy has its minimum at (f> = n), we see that when the temperature increases (/3 decreases), 
the total current becomes positive, and the minimum of the free energy shifts from (/> = tt to </> = 0. 

V. DISCUSSION 

We have studied in the previous sections the behavior of the Josephson current as a function of the spin coupling 
strength, and found that a tt junction behavior appears when this coupling is large enough. In view of an experimental 
realization, one must ask if the actual value of the spin coupling obtained with a given molecular magnet is large 
enough to observe this tt junction behavior. While a precise estimate, for a real molecule, of the magnetic coupling 
energy between the electronic spin and the molecular spin is beyond the scope of this paper, we can get a gross 
estimate by calculating the interaction energy of two magnetic dipoles at a distance typical of the molecular distance 
involved in our problem. Taking a spin S" = 10 for the molecule (as in Mnl2ac), and a distance '--^ 5 A, we find a 
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interaction energy O.lmeV, which is of the same order as the superconducting gap. This estimate shows that the 
TT junction regime due to spin coupHng may be reached experimentally. 

Let us now discuss some potential applications of our results. The system could beused as a Josephson current 
switch. Looking at the panel A of Fig. [71 we see that there is an abrupt change of the current sign as e goes through 
a specific value depending on the other parameters (it is e ~ 1.5 on the figure), while the current does not change 
much elsewhere. As e should be experimentally accessible using a gate voltage, a Josephson current switch could be 
implemented. Moreover, this implementation should be easier than in systems where the Josephson current changes 
sign several times as a parameter is varied. 

A more ambitious application would be to engineer a qubit with the system we describe in this work. Indeed, 
we have shown that, when varying some parameters, it is possible that the system behaves as a bistable junction 
(see for example panel B of Fig. [7]), where the system has a degenerate ground state. This feature can be effectively 
exploited to fashion a qubit system^, where the junction itself can be in a superposition^'' of the two ground states 
at either a phase difference of or tt. In contrast to the superconducting persistent current qubil^i, it is here the two 
phase states of the Josephson junction which provide the two states of the qubit. These qubits are therefore called 
superconducting phase qubits as in Ref.^3l- Similar to that in Ref.|4^, the coherent Rabi oscillations in our system 
could in principle be observed by a measurement of the phase sensitive sub-gap Andreev conductance across a high 
resistance tunnel contact between the qubit and a dirty metal wire^. 



VI. CONCLUSION 



To conclude, we have studied in this work the properties of the Josephson current between two superconductors 
through a single molecular magnet, which we modeled as a quantum dot plus a large frozen spin. We have shown 
that the coupling between the electronic spin on the dot and the molecular spin lead the system to behave as a tt 
junction. We have given a simple mechanism explaining this tt junction behavior, in terms ob bound state current 
and continuum current. 

We have shown moreover that the other parameters of the system give a precise control of this vr junction, allowing 
for example to reverse the tt shift and to bring the system to the normal junction state, or to an intermediate 
bistable state. This control of the tt shift can lead to useful applications, like a Josephson current switch, or could 
even be used to engineer a phase qubit. 

Possible topics of future study in such systems may include incorporating the dynamical nature of molecular spin— 
and quantum tunneling of the magnetization^, when the anisotropy barrier is not much larger than all the other 
energies of the problem. Also interesting would be to include electron-electron/electron-phonon interactions in the 
present work, in order to study the combination of the effects of the molecular spin and of the Coulomb charging 
energy, or/and the effect of electron-phonon interactions. 

Centre de Physique Theorique is UMR 6207 du CNRS, associated with Universite de la Mediterannee, Universite 
de Provence, and Universite de Toulon. 
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